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ABSTRACT 

This thesis contains some results of model calculations 
on the dynamics of point defects using a pseudomolecular 
model. Chapter I contains a discussion of the pseudomolec~ 
ular model in which a defect can be treated as a pseudo- 
molecule embedded in the rest of the crystal. The pseudo- 
molecular frequencies can be defined in terms of the effect- 
ive force constants obtainable in terms of the zero frequen- 
cy limit of the molecular Green functions. The self -entropy 
of a general point defect can te expressed as a function of the 
perturbed and unperturbed pseudomolecular frequencies by 
choosing the pseudomolecule such that the force constants 
between its boundary atoms and atoms of the unperturbed 
region remain unchanged. The express ions for self-entropy of 
substitutional defects and vacancies in some lattice models 
are given in Chapter III. 

In Chapter IV^a formalism for studying the time-depend- 
ent dynamics of point defects in terms of the spectral density 
matrices derivable from the molecular Green functions (for 
non-zero frequencies) has been given. The time-asymptotic 
form of a normal mode can be studied in terras of the isol- 
ated singularities and critical points of the spectral 
density functions. A sum rule for the spectral density 
matrices of the perfect pseudomolecul^ has been given. 



VI 


Neglecting relaxation of the medium around tlie vacancy, 
it has been shown that a vacancy in the simple cubic lattice 
does not have localized modes of vibration. A method of 
obtaining correlation functions of a point defect has been 
discussed. This method has been used to obtain an expres- 
sion for the Debye-Waller factor of cne of iiie nearest neighbours 
of a vacancy in the simple cubic lattice. 

For lattices with polarisable atoms a modified pseudo- 
molecular model in terms of the core and shell co-ordinates 
has to be used. Such a model has been discussed in Chapter 
V. It has been illustrated for a substitutional defect 
which does not change the' long-range Coulomb interactions 
in the simple cubic lattice. 



CHAPTER I 


INTRODUCTION 

The theory for investigation of the dynamioal proper- 
ties of imperfect crystals was given by Lifshitz ^ in 

connection with the optical properties of non-ideal cryst- 
als in the infra-red region, A similar method was later 
developed independently by Montroll and Potts Since at 
non-zero temperatures every substance contains a finite 
concentration of defects which play an important role in 
many physical properties like optical abso3rption, thermal 
conductivity/ diffusion etc./ the theoretical and experi- 
mental aspects of this problem were studied subsequently 
by many authors. The subject has been reviewed exhaustively 
in two recent articles by Maradudin^^^ . 

The earlier techniques used in this problem exploit 
the fact that a local perturbation on the masses and force 
constants due to a defect can be represented by a perturb- 
ation matrix having a dimension much smaller than that of full 
dynamical matrix of the perfect crystal. This enables a 
Green function method to be adopted for the evaluation of 
the perturbed eigenvectors and eigenfrequencies, in which 
the secular equation for the perturbed frequencies can be 
written in terms of a subspace/ the coordinate indices of 
which refer to those of the affected atoms. Such an approach is 
explicable Only to impurities which do not change the number of 

degrees of freedom associated with the perfect crystal. 

( 6 ) 

Stripp and Kirkwood' ^ have used perturbation theory for 
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A logical extension of the previous methods is to treat 

the defect region as a pseudomolecule which is embedded in 

the host crystal^ and which interacts with the latter. This 

approach is possible through a partitioning technique which 

enables one to write the exolicit equations of motion of the 

atoms in the defect region including its coupling to the rest 

(7 8 9 ^ 

of the crystal. This method has been used by many workers * ' ' 
in connection with different aspects of imperfect lattices. 

The effective interaction of the pseudomolecule with the 
host lattice can be evaluated exactly, for short range inter- 
actions, in terms of the Green functions and lattice parameters. 
In this thesis, the results of lattice dynamical studies of 
point defects on the basis of the pseudomolecular model will be 
presented. 

As an illustration of the use of this approach to the 
study of crystal defect problems mention may be made of 
computation of the vibrational self -entropy^ of a point 
defect. The self-entropy of a point defect at high temper- 
atures is expressible as a function of the ratio of the pert- 
urbed and unperturbed f requencies . As this problem 

reduces essentially to determining the perturbed spectrum, 

(13 14l 

the previous exact methods for determining the self -entropy' ‘ * 

are based on the Green function method or perturbation method. j 

( 11 ) ! 

Among the approximation methods, Huntington et al, have [ 

used the Einstein model for finding the perturbed frequencies 
of atoms in the vicinity of defect. They have calculated the 
contributions to the self-entropy, from the region of elastic I 
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( 15 ^ 

strain and surface effects separately. Theimer^ and Chandra 
etal'’ ' have assuned'Sie changes in freq;uencies due to the infest to be snail 
for finding the entropy of formation of a vacancy in the ionic 
crystals. Matthew^ has discussed the relative limitations 
of these approximation methods. Using the pseudomolecular 
model/ the self -entropy of a point defect can be obtained 
from the ratio of the products of frequencies of the perfect 
and imperfect pseudomolecules . This method is physically 
more realistic as the exact interaction of impurity with the 
host lattice is taken into account. Furthermore, the pseudo- 
molecular method has the advantage of being applicable to 
a general point defect which may involve a change in the 
number of degrees of freedom of the lattice. Computationally/ 
also, this method is not more involved than the conventional 
methods as the secular determinants of the pseudomole- 

cules can be blpck-diagonalised by using the point symmetry 
of the lattice. In the present model calculations we have 
obtained the pseudomolecular frequencies for substitutional 
impurities and vacancies in the square, simple cubic and 
B.C.C. lattices. 

The pseudomolecular method is also useful for the study 
of detailed dynamics of point defects, a knowledge of which 
is required for studying physical properties like optical 
absorption and Mossbauer line-width etc,^"^'^^,’ This method 
enables the exact dynamics of impurity to be studied in terns 
of the spectral density functions of the normal modes of 
vibration, of the defect pseudoraolecule— The spectral density 
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functions are qiaite complicated, in general, and therefore 
the integrals required for obtaining the explicit time- 
de-'endence of the atomic displacements are difficult to 
evaluate even for Isotope defects in the linear lattices, 
Hov/ever, a knov/ledge of the time-asymptotic -behaviour of the 
displacements, which is governed by the nature of spectral 
density functions at the critical points and singularities, 
can be used to investigate the possibilities of occurrence of 
persistent modes of vibration. Relaxing modes of vibration 
correspond to the well-behaved part and persistent modes to 
the isolated singularities of the spectral density functions. 
The present method shows exactly that a vacancy in the simple 
cubic lattice (relaxation around the vacancy has not been 
taken into account) does not have localized modes of vibration, 
a result which was proved earlier by Lengeler and Ludwig^ 
by an approximation method, A formalism for computing the 
correlation functions of a general point defect by the pseudo- 
molecular model has also been given. 

For a lattice with polarisable atoms. Page and Strauch^^^^ 
have given a perturbation method based on Lifshitz's Green 
functioA. method for studying the dynamics of an impurity 
which does not change the long-range Coulomb terms. A pseudo- 
molecular model for such a case has been discussed in 
Chapter V of the present work. 
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CHAPTER II 

A PSEUDOMOLECULAR MODEL 
2.1 Par-titloninq Technique 

The dynamical equations of atoms belonging to a certain 
subspace of a crystal can be separated by partitioning various 
matrices occurring in the equations of motion of the entire 
crystal. The procedure is similar to what has been used by 
Lowdin^^*^^ for contraction of secular equations in connection 
with quantum-chemical problems. Lowdin has shown that a sec- 
ular equation in a certain representation can be contracted 
to a similar secular equation of smaller degree and referring 
to a subspace of the original basis by partitioning the matr- 
ices suitably. The contracted secular equation so obtained 
contains energy as a non-linear parameter. Inspite of the 
simple nature of the original secular equation, solution of the 
contracted secular equation by iteration procedures is more 
accurate, as has been pointed out by LlSwdin^^^^. Also, the 
solution of the contracted secular equation by Newton- 
Raphson*s second order iteration method assumes a very simple 
form which shows its analogy with the variation method and 
connection with the general perturbation technique. 

The equation of motion of a perfect lattice having N 
degrees of freedom can be written in matrix notation as, 

= L X {t) 


M X (t) 




( 2 . 1 ) 
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A A . 

where L is a (NXN) matrix of atomic force constants, M is 
a H-dimensional diagonal matrix of atomic masses and X(t) 
is a N-dimensional vector having time-dependent atomic 
displacements as its components. 

Assuming X(t) = U(^) in equation (2.1), the 

following time-independent equation of motion of the latt- 
ice is obtained. 


(L-MO^) ) = 0 


( 2 . 2 ) 


We partition the matrices L, M and vector U(<^ ) in eqn, 
(2.2) as follows. 


A 

M 





(2.3) 


where matrices C, m and vector 22^^'^ refer to a region 
Q of the perfect crystal which will subsequently be 
perturbed by a point defect. The corresponding matrices 

A A 

for the imperfect crystal shall be denoted by L’ , M‘ and 
U’, respectively. The partitioned forms of these perturbed 
matrices shall be similar to those given by eqns.(2.3) with 
all the submatrices, except A and ?5 ;a/ as primed. and 

2 ^(co) refer to the unperturbed part crystal. ^ 

and are matrices of the force constants between the 

and Q. We assume that the region 


a ton® of regions 



has n and n' degrees of freedom associated with it for 
the perfect and imperfect crystals, respectively. 
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Using eqns. (2.3) in eqn. (2.2) we get the following two 
sets of N~n and n equations in matrix notation. 


(A-MjSj^ta2)Uj^+BU2 = 0 

(2.4) 

B%^+(C-m = Q 

(2.5) 

Prom eqns. (2,4) and (2,5), the equations for and 

referring to the regions P and p , respectively, 

1 ^ 

, be separated as. 

(i-%ca 2^1 6^2) = 0 

(2.6) 


(2.7) 


2,2 The Pseudomolecular Model for a Point Defect 

A pseudomolecular model for a point defect can be 
given from the contracted equation (2,7). Eqn. (2,7) can 
be interpreted as the equation of motion of a pseudomole- 
cule having n degrees of freedom and embedded in the rest 
of the crystal having (N-n) degrees of freedom. Comparing 
eqn. (2.7) with the original equation of motion viz,, 
eqn. (2,2), the additional term B'^(^-I^ji^ooi^)§ here can be 
asguraedl to re^re^eht an 0f£^iitiv4 iriteraction oi the • ifSeiado* 
molecule with the rest of the lattice. This type of appro- 
ach has been used by many authors for various lattice- 

dynamical problems. 
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Litzman and Rozsa^"^^ have calculated the frequencies 
of localized modes of a defect in the square lattice by 
taking a pseudomolecule of thirteen atoms. They have part- 
itioned the total space into three subspaces. Explicit 
evaluation of ”^B matrix, as has been mentioned 

by Litzman and Rozsa, is an extremely difficult task except 
for the case of a linear lattice. However, by using certain 
identities satisfied by the submatrices of the matrix 
(L-M tu2) and its inverse it can be written in 

terms of the known matrices. 

AA, O /'-/s 0*****i 

By partitioning the matrices and ) 

as shown below 


(L-Mu2) 


bT 


A 

B 


C-m o' 


\ 


( 2 . 8 ) 


2-1 f ) 

^ 0^(0 ) 

,(7) 


Ge(<^ ) 
Gci^ ) 


it can be proved 

G~ ^ ( oJ ) = C-m «2 _bT ( O 2 ) “ ^ 


(2.9) 


Substituting eqn. (2.9) in eqn. (2.7), the dynamical 
equation of the pseudomolecule can be written as 

=0 ( 2 . 10 ) 
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By using the notation explained for eqn. (2.3)^ the 
dynamical equation for the imperfect pseudomolecule can 
similarly be written by priming the various submatrices 
occurring in eqn. (2,7) as, 

(C' -m' ) )u^(« ) =0. (2.11) 

To evaluate B''^(A-M^ o:*^) ^B' matrix in eqn. (2,11)/ 
we use the following property of this matrix. A typical 
element (t,?*) of the matrix b'^(A-M^ oJ^) B is given by. 


B^(A-Maai 2)-1b = ^ bT(^ ,R) (A-M^ £>^2) B(R‘4'). 

Jh R.R» 


( 2 , 12 ) 


In eqn. (2.12), and 4' refer to the indices of the pseudo- 
molecule and ,R and R' refer to those of the unperturbed 
subspace i If the force constant between an -iatom and a 
R atom is denoted by oqn. (2.12) can be written as 





ill All the R All the 

atoms inter- R’ atoms 

acting with interact- 

i. ing with 4' 




( 2 . 13 ) 


It follows from eqn, (2,13) that for defects which 
do not involve a change in the degrees of freedom asso- 
ciated with the crystal/ B*'^(A-M;i^w2) ^b' and b’^(A-Mjii^o3 2)B 

I 

will be exactly identical if 'T^p,,s do not change, For 
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short range interactions, this can always be done by 
choosing the pseudomolecule in such a \fay that the subst- 
ituted or removed atoms do not directly interact with the 
atoms of the unperturbed region. But if the defect changes 
the number cf degrees cffreedomof the crystal (i.e., n' / n) , 

Ov— 1*^ . Am ^ O . ^ 

B* and B^(A-M^c^'"^)B shall differ in 
dimension , However, with the above-mentioned choice of 
the pseudomolecule, these matrices will have exactly ident- 
ical non-zero elements. Therefore B ’ '^(A-M^ib2) ' can be 

obtained from b'^CA-Mco^) ^3 addition (for n'> n) or by 
deletion in' <. n) of rows and columns of zeros correspond- 
ing to the indices of added or removed atoms , respectively. 


B'^(A-Mj^«2)“1b' = ^ ,(2.15) 

where 1 denotes contraction (or addition) for n'< n 

L -JC/a- 

(or n'>n), respectively. 

It is possible to derive some useful relations between 
the secular determinants of the perfect and imperfect 
crystals . 

The secular equation of the perfect crystal can be 
written from eqn. (2.2) as. 


L-M CO 



0 


(2.16) 


By partitioning the matrices in eqn. (2.16) as shown by 
eqn. (2.8) we get 
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A 2 
L-McJ 


0 , 


,A A 2 

A 

A— 

b 

bT 

C-m<.2 



. bT 

C-mi-J^I 


0 \ 

bT 

A A 2 

C-m''-’' 

As /V ry 

|C-mo:> j 


A 

0 B 


, 0 0 
A A 


A-M; 


o2 


I+l 


.V A A , — 1 " m A 9 


-1 


A 

'O B 


iO o- 


0 


(2,17) 


Similarly, the secular equation for the imperfect crystal 
can be written as 

j L ‘ -M ’ co2 I =1 )6 ‘ -m’ 1 1- (C ’ "^b ' ^ , j 

(2.18) 


Prom eqns. (2.17) and (2.18) we get 

!£t,^i.,(o2| ^ ^ 2 ^ -ll , 1 

!C-.rhco2-BT(i4)i^co2)-lB | * ^2.19) 


Eqn. (2,19) shows that the ratio of the secular 
determinants of the imperfect and perfect crystals is 
given by the ratio of the contracted secular determinants 
of the imperfect and perfect pseudomolecules. This result 
will be used in Chapter III to find the self -entropy of 
point defects. 
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9 ^ ^ Effective Interaction Matrix 

The pseudomolecular method reduces the degree of 
the secular equation and facilitates inclusion of any changes in 
the number of degrees of freedom associated with the perfect 
region due to the impurity and it tabes into account the 
exact interaction of the defect with the crystal. Equation 
of motion of the atoms belonging to regions ( ^ and | 2 
have been written down separately as eqns. (2,6) and (2,7), 
respectively. These equations are exact as the interaction', 
between the atoms of ^2 'baleen into account by 

" . /s 2 \ “"^oT o 

the effective interaction matrices B (.C-m ^ ) b ana 

Some interesting results can be derived from the 
effective interaction matrices. By using approximate 
expansions of (A— M^o^) , ICrurrihansl and Matthew have 

obtained equations for determining the frequencies of very 
high-frequency localized modes and low-frequency resonances 
from the secular determinant of the impurity pseudomolecule. 
These equations have been claimed to be as accurate as some 
of the other numerically computed results based on some 
special choice of the force constants. 

In the present worb^we have not considered the scat- 
tering of lattice waves by general point defects. If the 
"^ubspace U 3 ^.C^ 2 , 1 ) is defined such that it refers to the 
,^ices of unperturbed atoms and the three centre of mass 
^dinates of an atom having the mass of the molecular 
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impurity, the resulting equation (2,5) will be the equation 
for a lat'^’ice having the same mniber c£ d^rees of freedom as the 
perfect crystal and can then be solved by Lifshitc’s method. 
Using this method, Wagner^ has studied the scattering of • 
phonons by molecular impurities. By eliminating the addi- 
tional degrees of freedom due to the molecule, he has sho^-m 
that the usual local perturbation matrix of the imperfect 
crystal contains an ad'^itional tenn having a form similar 
to B(C-m«^)""^ b”^ which has poles at the molecular frequen- 
cies; (C-moo^)”^ being the molecular Green function matrix. 
Wagner has further shown that the t-matrix of the scatter- 
ing formalism and, conseqaently^ the scattering cross-section 
show sharp resonances near the molecular frequencies, A 

similar aoproach regarding the defect as a particle ertibed- 

( 22 ) 

ded in the phonon field^has been given by KcuiThansl 

for studying the scattering of lattice vibrations by vacancy- 


type defects. 

Arn,/\ ^ O. — 1 / ;? ^.«2^~l AT 

The matrices B^(A-M^cO'“) b and B(-^-m£^ ) B can 
be computed in terms of the Green function submatrices G^ 

^ 1 A — 1 

and G^ from eqn, (2.9). G,^ is not difficult to obtain 

as it has the point symmetry of the lattice and can there- 
fore be block-diagonalised by using group-theoretical methods. 
In cases where group-theoretical methods of inverting matri- 
ces by block-diagonalisation become involved an alternative 
method suggested by Litzman^'^^ may be used. This method 
is based on a theorem which states that a unitary matrix 
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constructed from the kno^\m eigen vectors of any matrix 
partially diagonalises all the matrices commuting with 
that matrix. 

The pseudomolecular model presented in this Chapter 
has to be modified for lattices with long-range interactions 
due to polarizable atoms. This will be discus.-ed in detail 
in Chapter V. 
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CHAPTER III 


VIBRATIONAL SELF -ENTROPY OP POINT DEFECTS 
3.1 I n tr oduc tion 

The vibrational self -entropy of a point defect As, 
defined as the thermal entropy of formation of the defect, 
in the harmonic approximation is given by^^^^. 


= S'-S 


v 

r 




1-exp (- +(!l£rN 

kgTj T ^ 


[_®xp( I 


(3.1) 


i 


i-exp(-^i)l r exp(e^)- 




j 


where S‘,oj' and S,co) are the entropy and frequencies of the 
imperfect and perfect crystals, respectively, kg is the 
Boltemann's constant and T is the temperature of the crystal. 
For high temperatures (i.e., kg Tyy'h^j^), eqn. (3.1) 
reduces to 


AS ^ 


I 




-N) 





+ (N‘-W)k3, 

(3.2) 


where N’ and N are degrees of freedom associated with the 
imperfect and perfect crystals, respectively. In eqn, (3.2) 
terms of order > ( £, ...- ^ ) have been neglected, being 

^ JCg J- ^ 

the Debye frequency. 
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From eqn. (3.2) it is clear that knowledge of the 
frequency-spectrum of the perturbed and unperturbed 
crystals is required for evaluation of JHaSi This problem 
has been studied in detail by many authors ^ . A series 
of exact calculations of A S for a mass defect in the 
simple cubic lattice/ a vacancy in the f.c.c. lattice, 

substitutional defects in simple cubic lattices have been 

. ( 13) 

given, respectively, by Naradelli and Tettaman'^i , 

Nardelli and Terzi^^'^^ and Mahanty^^^^, These calculations 

are all based on Lifshitz’s Green functions method. It can 

( 23 ^ 

be shown that for N‘=N, AS is expressible as , 

AS = - limit M|»A(«) (3,3) 

( 0—^0 2 

... A A 

where A(o5^) == [ I+G P j with G and P as the Green function 
and perturbation matrices of the perfect and imperfect 
crystals, respectively, is called the defect determinant. 

A 

The perturbation matrix P has only a few non-vanishing 
elements which can be grouped together to form a submatrix 
p in a certain subspace r 2 of the crystal. Furthermore, 
the submatrices p and g have the point symmetry of the sub- 
space feature can be exploited to block-diago- 

nalise A(co) by group-theoretical methodswhich enable 
one to determine the contribution*, to AS ,from the modes 
belonging to various irreducible representations of the 
point symmetry group of ^ 2 * have given the block-^ 
diagonalised form of A(m) for a substitutional impurity 



17 


in the simple cubic lattice in Appendix I. Benedek and 
Nardelli ^ have expressed in terms of the numerator 
and resonance-denominator (= (oj ) ) of the scattering 
matrix. 

One of the earliest approximation methods of deter- 
mining ^ S of a substitutional defect based on the pert- 
urbation theory was given by Stripp and Kirkwood^ ^ . 

Another more commonly used approximation method based on 
the Einstein approximation (regarding each atom in the 
neighbourhood of defect to be oscillating in an average 
field produced by the rest of the atoms which are held 
stationary) has been given by Huntington et al.^^^^. They 
have obtained S of a vacancy in copper by using the 
Einstein approximation for calculating the contribution of 
imm.ediate neighbourhood of the vacancy to Zl S and^ elasticity 

theory for that of the rest of the region, Nardelli and 
( 2 ^) 

Chiarotti have calculated the entropy of a vacancy in 
crystals of inert gases by using a similar approximation 
method. For alkali halides, Theimer^^^^ and Chandra et 
have expressed £lkS in terms of the change in trace of the 
dynamical matrix of the perfect crystal due to the defect 
assuming the changes in frequencies due to the defect 
' /^e small). Matthew^ has discussed the relative reli^ 

' .^ty of these approximation methods by comparing the 
‘ > ts from them with the exact results for substitutional 



in linear and cubic lattices. 
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Recently^K. Mukher jee ^ has given a quasi- 

thermodynamic method of obtaining ;;^S of vacancies in 
f.c.c. metals. This method is based on the exoerimental 
fact that the equilibrium concentration of vacancies in 

_4 

f.c.c. metals near the melting point lies between lo 

— ( 27 ) 
and lo Subsequently, K. Mukherjee has applied 

this method to solid argon and the results obtained are 

in good agreement with experimental values within the 

approximation of the Einstein model of solids. 

We shall discuss the pseudomolecular model of finding 

the self-entropy of a general point defect^^^^ in^3.2 and 

^3.3^and some model calculations based on it in ^.4, ^3. 5 

contains a discussion of the various theoretical and 

experimental methods of finding As. 


^3,2 Method Based on Pseudomolecular Frequencies 

An exact method of finding the self -entropy of a 
general point defect based on the pseudomolecular model 
|discussed in Chapter Il)has been given by Mahanty^^*^^ , The 
expression for As follows from the use of eqn, (2.19) in 
eqn. (3.2). 

Using the notation of ^2.1 in eqn. (3.2), can 


be written as 
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AS = - 'Ll IM 

-T-^^ 'l^T' ,i,l iflV 


+ (N*-N)k, 


■B 


T 


In 


(— ) 
IcbT^ 


IC'-B'’’ 

B'i 

1 rri 1 

[S-bT ^-1 

B 1 

im'l 


_ limit 


-^L 

2 


t (]SI'-N)kB 


, n |c--§tT A~^ Im' ! 

kB^' |-a_%T l-lg _ flij2 1 |5i.| I 

n' „ J_ 


ti ^2(n*~n) I 

Krr^/ — 


ff, 


,-f(N‘-N)kg 


m i 


'1%*' 


-5 ( 

TT 

^=1 


m' 


+ (N‘-.N)kB 


(3.4) 


liere cO^ are given, respectively, by the roots 

€ 

f the following two equations 


C' - B-’^ A -1 B' _ fi'ca^ 


= 0 


(3.5) 


id 


C -• A B ^ m 


limit G”^ (w ) -m o3 2 


= 0, 


cO ^ 0 

lere use has been made of eqn, (2,9) in eqn. (3,6), 

Equations (3,5) and (3.6) give n* and n frequencies 
and respectively) associated with region ^2 
le imperfect and perfect crystals. These frequencies may 
! called the pseudomolecular frequencies as they have been 


(3. 6) 


■itten in terms of effective force constant matrices 

B* and b'^ B (of dimensions n^ and n , respectively) 
ich are obtainable from the effective interaction matrices 
the imperfect and perfect ‘pseudomolecules (Chapter II) 
the zero frequency limit. The self-^entropy of a point 
feet is jbheref ore /expressible as a function of the perturbed 
i unperturbed pseudoitKolecnalar frequencies in Idie sense 
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'^3.3 Divergence Problem in Determining Gc 

The matrix elements of ) occurring in eqn. (3,6) 

are defined by the following equation^ 



k 

where (k) is the component fof the k^^ eigen vector of 
the dynamical matrix and ,:.o^(k) is the corresponding eigen- 
frequency, The summation in eqn, (3.7) is over all the 
vibrational modes of the perfect lattice. As mentioned 
earlier, G^(o ) has the point symmetry of the crystal and 
can^ therefore, be block-diagonalised by using group-theoret- 
ical methods and then inverted. A divergence difficulty 
arises in this method for the limit cj—*o. From eqn. (3,7) 
it is obvious that the roots cd2(k) = o will give divergent 
elements of Green functions for (O — ?0, Such roots, which 
correspond to the translational modes of the lattice, are 
always present in the secular determinants of perfect and 
imperfect crystals. These 53 ^(k) = o roots persist even in 
the contracted secular determinants of perfect and imperfect 
pseudomolecules . 

To overcome this difficulty we construct a unitary 
matrix V which contains vectors corresponding to the trans- 
lational modes, and other orthogonal vectors such that 

, A. Afc 

Gq V is block-diagonalised. Such a V matrix can always 
be constructed by taking linear coiDbinations of the syirenetry- 
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adapted vectors belonging to various irreducible represent- 
ations (IR) of the point-symmetry group of subspace T 2« 

It is found that with this particular choice, matrix 

/'+ A . A 

V )V contains a divergent block corresponding to the 

translational modes. Other blocks, corresponding to the 
vectors orthogonal to the translational modes, involve such 
differences of the elements of that divergences in 

them cancel out. The divergent terms can generally be 
simplified considerably by using the sum-rules satisfied 
by Green functions ^ in the zero frequency limit. 

^3,4 Calculation of the self-entropy of a substitutional 
Defect and a Vacancy in Some Lattice Models 

We shall give here sane results of computations of the 
pseudomolecular frequencies of substutional impurities and 
vacancies in the square, simple cubic and b.c.c. lattices. 

The interatomic interactions in the lattices are assumed to 
be of the nearest-neighbour type. Relaxation of the structure 
around the impurities can be include;^ in the parameters of 
the defects in the model calculations given in ^3.4a-^ 3,4c. 

^3.4a. The Isotropic Square Lattice 
(i) Substitutional impurity : 

This hypothetical model is considered here to illustrate 
the details of a typical calculation. For mathematical 
simplicity each atom is assTomed to have only one component 
of displacement so that polarisation effects can be ignored; 
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A substitutional impurity having mass m' and force 
constant t‘ is taken to be at the origin of subspace 
r 2 (fig (lb)), m and t are the ness and force constant, respect- 
ively, of a normal atom (fig (la)). Here n' = n = 5, 

The integral form of elements of Gq(<*®) is given by^^^^» 

g(i.,i';<^) = g(%-4/ 


= 1 f f Cos(( il“-^)ki) Cos((i2-^2)lc2) dKldks ^ 3 ^ 3 ^ 

( 21 T)^ Q Q 2 IT ( 2 -cos 3^ -cos k 2 )-mcu^ 

In this isotropic model, g(|., |.‘;co) = g (4/ -4 7 ^). The 
indices needed to specify the elements of are 0,1, 2,3,4, 
where O refers to the impurity site and 1, 2,3,4 are the 
nearest neighbours taken in a cyclic order, this group of 
atoms constituting the pseudomolecule. 

The Gq{<-'3) matrix is given in Appendix II. The matrix 
of symmetry adapted vectors i of the point symmetry group 

A 

D 4 of f 2 * matrix V constructed out of the vectors of 

§ as explained in ^3.3 are given in Appendix III. 

The Sum-rules satisfied by the Green functions are. 




Isotropic square lattice 



Fig* la Pseudomolecule for the perfect crystal 



Fig.Ib Pseudomotecul(^ for the square lattice cohtaNngi 
a substitutional impurity at the origin. 
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The block-diagonalised. form of Is 


-+5g(l,0;«) 


V^G^(tXi)v= 


O ^+(g(2,0;«) 

2 ^ -g(1.0;c^)) 


(g(o^o?w) 0 

-gU^O;^^ ) 


o (gCo#o,*‘^) 
-g(2,0;c4) 


^ (3.11) 

In the limit w — ^ 0, the divergent term occurs only 
in the first block which corresponds to the translational 

A 

mode written in the first column of V. 

Inverting eqn, (3.11) and performing inverse trans formation ^ 

A-1 

the G(-. matrix obtained is, 
limit (^ )E C-B^ A B 


-r(n2,6 4-4) %i 2-10 
2(4-T!)(n-2) 2(4-»0 


Y( n2..2rt-4) 
2(4-Ti) (lT -2) 


Y(2-r0 

olA-si') 


2(4-n) 


~^(Ti2~6hHi-4) 2~T0 

2C4-IT) (ti-2) 2(4-1T) 


'!rf'H^~6»"-4) 
2(4-T0 (tT-2) 


YfT>2~2If-4) 
2C4-TT) (■n~2) 

•^(24?) 


^Y(r?2..6n>»4) 


2(4-^¥) 


-y(Tr2^2ff»4) V€2~T0 


t-K} CP- 


T(2-TT) 

2(4--n) 


where use has been made of the sum-rules given by 
eqns. (3,9) and (3,lo) in the limit co — ^0. 

The sum of^each row of in eqn. (3.12) is zero 

showing thereby that the pseudomolecule has a zero frequen- 
cy associated with its translational mode. 

The matrix C of force constants between atoms of the perfect 
region is^ 


i 4 3' -T 

/ -It 4'5- 


/V 

C 


-H 


o 



\ 


-T 




0 

0 


0 43' 

0 0 



(3.13) 


Matrix B obtained by subtracting eqn, (3.12) from 

eqn, (3.13) is^ 


A B= 


O 


0 4T + 


O 


0 


t (T(2^6i?+4) ( u-2) 

2(4-T!) (Tr-2) 2(4-Tr) 


-Y( fi^-2Tr-4) -2f(Tr~2) 

2(4-10 (r-2) 2(4-7T) 


'^(n‘-2) .'.,"y(n^-6<r t4) 1?(TT-.2) -'T( tT^-2<T-4) 

2(4- ) ^ ^ 2(4-u)(r-2) 2(4-?T) 2 (4-Ti) (T?--2) 

T( P^-2iT- 4) h(TF-2) ^(n2_6vf4^4) ^(0-2) 

2(4-fr) (ir-2) 2(4-n) 2(4-7r) (n-2) 2(4-it) 


ir(iT-2) 

2(4-11) 


-l(T<^-2<T-4) '<'(tr-2) ,i 

2(4-T0 (5T-2) 2(4-TT) 2(4-7f) (mW 


(3.14) 
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The matrix in eqn. (3.14) has a row rind a of sero etensite 

corresDonding to the index of interior atom o which does not 
interact with any of the atoms of region 1"’;]^. This follows 
from eqn, (2.13) for the elements of B A B in the nearest 
neighbours interaction approximation, 

A ^ 

The mat'^ix C' of force constants of atoms belonging to region 
f 2 pertttrbcd by a substitutional defect^ is given as# 


C’ = 


4T’ 


-T’ 

-X' 

-y' 

3Y+T' 

0 

0 

-Y' 

0 

ST+q* 

0 

-V 

0 

0 

3r+Y' 

-Y* 

0 

0 

0 



(3,15) 


In this case matrix b'^‘ "a*”^ B* is exactly identical with 
the matrix a“^ B given by eqn. (3.14). 

The perturbed pseudomolecular frequencies obtained by 
substituting eqns. (3,14) and (3,15) in eqn. (3.5) are/ 


. I _ o. 2 _ 4m4-m* . 

CO 0/ 0^ I 0 / 

o 1 mm 


,2 _ q* ^ 2y(H-2) 

^2 " ^ m (4-Tr) 


and 




■2 


m m(!r-2) 


(3.16) 


The unperturbed psueoraolecular frequencies can be obtained 
by solving eqn. (3,6) or by taking the limit f* — fS and 
m‘-~> m in eqn, (3.16). These, therefore, are given by. 
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and 



TT ' 


-2 ^ TIT 
'"*'2 m(4-7T) 


2 _ , 2 _ TTil 

*^"^3 ^ 4 m(u-2) 


(3.17) 


Substituting eqns , (3.17) and (3.16) (dropping out 

, ~ eqn. (3.4) the self entropy can be 

"^o o 

written down in terms of parameters T ' ^ T /in' and m of 
the lattice. 


( ii) Vacancy in the Square Lattice 

In the scalar model considered here a vacancy decreeses 
the rurriber caf (3egrees of freedom associated with by one (i.e.^n 
-n’ = 1) . Matrices C' and B’ are therefore (4x4) 

matrices. If the vacant site is taken to be at the origin 
of region f 2 / matrix B' a” B' can be obtained from 
gT ^-1 g (i.e., eqn. 3.14) by deleting the first row and 
first column of zero elements (^2.2). 

A , 

Matrix C‘ is given as, 

C' = 3TfI, (3.18) 

where 

I = (4x4) unit matrix. 

The pseudomolecular frequencies in this case can be 

. ^ , (10) 

shown to be , 


2 


== O; 


2 _ 2? ,5r-2^ 

^ A- m \ jr ‘ mf 

2 


2 2 27: 

^ ~ ^ “ m(Tr-2) 


( 3 . 19 ) 
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3 . 4b The Isotropic Simple Cubic Lattice 
(i) Perfect. Lattice ; 

Consider a simple cubic lattice consisting of atoms 
of mass m each and nearest neighbour force constants T . 
The atoms are located at the sites given by 

= i^ai -ii- (3,20) 


v7here 


^2' ^ 3 ^ ^ lattice constant. 

An atom at the origin (O/O/O) and its six nearest 

neighbours at lattice sites (l/ 0 ,o); (o,l/ 0 ); ( 0 , 0 / 1 ); 

(1,0,0); (O/l/O); and (0,0,1), as shown in Fig. (2), shall 

be taken as the atoms constituting a pseudomolecule. The 

elements of Green function* matrix G(^(co ) are given by the 

( 28) 

following integrals ^ ^ , 


g(l,i«.co) = 


2rr 2ir 2^ 

( 


(2rt )3 i j 

000 





2T(3-Cos k^-Cos k 2 “Cos k 2 )-mco 

(3.21) 


Using group-theoretical methods, the subspace \2 can be 
reduced into the following IRs of group 




(3.22) 


V and 6 matrices are given in Appendix IV. Appendix IV 
also contains and the matrix's’^ A obtained from 

them. Using these in eqn. (3.6), the unperturbed pseudo- 
molecular frequencies obtained are. 
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Fig . 2 Pseudomoiecule for the perfect isotropic 
simple cubic lattice 3*4b)» 
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c.\ = 


= 




0- 

u, c>>2 RT 


^<^2 2 _ 

CO. = (.o = — 

4 5 m 


2 

&» = _ 

7 m 


A, - IR , 
Ig 


( 3 . 23 ) 


rg(o,o,o^g("2';o;d)T ra'^lu-™' 


( 3 . 24 ) 


j^g( 0 ,o,o) +g( 2,0/0) -g(l/l,o)] Sg-IR, 

( 3 . 25 ) 


where the numerical values of Green functions computed 

/ ^ w \ 

by Fourier-series method' have been used. 

ii) Substitutional Defect at the Origin ; 

For a substitutional defect of mass m* and force 


constant t’ placed at the origin of subspace 

C’ has the following form 

l\ 

i . 

n 

' 2 ' 

\ 


7’S* 


-Tf’ 

-y 


-n' 



1 -IS' 

1 

6Y+-’S' 

0 

0 

0 

0 

0 ' 


1 

0 


0 

0 

0 

0 

■’N. 

c = 

-V 

0 

0 

6 Ti tV* 

0 

0 

0 


-r 

0 

0 

0 


0 

0 


i 

0 

0 

0 

0 

61 

° 1 


\ 

0 

0 

0 

0 

0 

63 

( 3 . 26 ) 

The matrix 

A i A 

A“^ B* 

is identical with 

matrix 

B-^ A 

B 


given in Appendix 17 , 



The perturbed pseudomolecular frequencies are 


' ^ /2 _ 6m+m’ 

- 0/ <.o„ j— 1’ 

1 2 rnm' 




'2 _,,.2 _ 1 
3 4 m 


g{o, 0 , 0 ) -g{ 27 o,o) 


T +3.762^ 


m 


P, -IR, 
lu ' 




i2 ^ ,2 


Ca^ 


m 


. ^ ^ — X f 

g(0,0,0)-tg(2,0,0) -2g(l,l,o) J 


-(-4.417 ^ 

m 



(3.27) 


(iii) Vacancy at the Origin ; 

For a vacancy at the origin (0/0,0) , matrix C' can 
be written as 

c' = STSi, (3.28) 

where I is a (6x6) unit matrix. 

/N A ^ t . 

B' A B ‘ IS a (6x6) matrix and is obtained from 
b"^ a ^ B (Appendix IV) by deletion of the first row and 
fiiCfft column. 

The eigen-frequencies obtained are. 


2 

= o 




II 2 = Q 


..2 


CO 


• I 2 


• > 


_ 1 i 2 _ j. r_ i 1(1 

4 rn Lg(o7o7o)-g(27o7o) J 


- 3,762 


m 


: F- -IR, 
lu ' 


6 


= 4 - 

m t' 


g(o,o,o)+g(2,0/0) -2g(l 


T 

/liO) 


(3.29) 
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It may be mentioned here that each set of frequencies given 
by eqns. ( 3. 26) - (3, 28) occurs three-times as the X, Y and 
S comnonents of atomic displacements are decoupled in the 
scalar model considered here. 


3,4c The B.C.C. Lattice 

The model chosen here is the one that has been used 
earlier in connection with the phonon scattering^ . The 
nearest neighbour central and non-central force constants 
are and T 2 * The equations of motion are written in 
terms of the constants and &2 which are defined as. 


h 


+ r 


fi = 




(3.30) 


V matrix used here is given in Appendix V. It contains 
three transitional vectors corresponding to translations 
in X, Y and Z directions. 

The calculations have been done for a substitutional 
impurity. In the limit ^2 results have also 

been obtained for a vacancy. In this limit the three 
components of the displacement of each atom get decoupled 
in a manner analogous to the previous example of the simple 
cubic lattice. 

The unperturbed pseudomolecular frequencies are given 
by the roots of the following equation. 



33 




(-mCo2)x lOT 




t»j 2 




/Aig 


/l _ mco^ '{ ^ X ) m^-ca^-m (. ^^ i ^ . - 2 )+ — x/ ^ -mco^) 

\P /Eg I KIvI~2l 2 :. KM-2L2]2P2g^\^i / 


( 


( 1 


X — -m (W 


2\3 


„ 1 X / ” -mw^ ) 

^*' 2 u ^ '^''^9 




(3.31) 


are given 


ad the perturbed pseudomolecular frequencies 
7 the roots of the following equation: 

iL "'(e 

lu ^ 

I KM-2L '1 I ^ m-2L^ 


( Aji^- Aji^) A Pj -2 ^ (I 


/a, 


lu 


i^af+A^-mw^V =0 

•'hi 


ig 


(3,32) 
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‘‘ 21 JU 


VS , 2 

— Am<h> 

9 


. S ' 

55^72 ^ 3^2-3r2Af; 

- J ^^2 - 3J2 A^2 


(3.33) 
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with Am = m'-m? |3^; = |-2 - |^2 

( 31) 

the notation of Green functions used earlier / 

^ ^ ^5Do“^D4^'^^5D32“5D2l)‘*2(gD22~5D31^'^2(gN4'^gN3) ; 

F = ^5do“9d4^ +^5d32“5d21^"‘‘^'?D22“9D31^ “9N3"5n4' ?' 

^ = ^5Do“5D4)'‘"^gD21“?rD32^''‘^gN3^9N4) * 

^ ~ ^*^Do“5d4^‘*'^^D32“^D21^ '^^^^D22'"5D31^ ' 

^ “ ^^Do"^D 22^ "^^'^D4”'?D22)'^^5D21"5D3l^'^’-5D32"9D3l)‘^2(g^3-gj^4) / 

^ ~ ^^Do“%21^'^^5d4“^D32^ "%4'‘%3 ; 

^ = ^gDo-'5D4)-^^gD21-^D32^"%4-%3; 

^ = ^%o"5d22^'^^^D21-^D22^'*‘^5d4-9D31^^^5D32"^D31^'^^^4“%3^' 

u = ■! ■|^9(gQQ-gDi)+(gQ4~gjj^) + (gp2i”9Dl)‘^2(gj322“gDl^‘‘“^^^D31~%l^ 

'*‘^^32“9di)^ ; 

= ^%o'-5d21^'*^^%4“5D32>-^%4-%3^' 

= f ^ ®9 ni“3n4-9n3 ) • 

In the limit ^2 "'^0 and ^1 ="!r / all the non- 

diagonal Green functions g^jj^/ gjj 3 and become zero 
and gD 2 i = gD 33 _ = gD32^^^^* functions 

satisfy the following sum-rules, 



35 


S^Do-S^Dl " 


®gDl-^gDo‘^D4'^^^D2l'^^^D31^ = 0 ^3* 35) 

Using cqns, (3.35) in eqn, (3,32) the secular equation for 
the perturbed pseudomolecule becomes/ 

(_w 2^3 x5(in4A§)j 9T‘-(mt§Am) ft>2 I + |_( A m) ^ <^2 

I. 9 i_ 9 j 81 j 


Si / 1 9?^ \ 1 2( 

xj~+T-r-m//l X + T'-lf -mc.o^'l X i 1 -|"T'-T-m CO 


X 


= O 


r J J 

(3.36) 

The secular equation for a vacancy in this model is given 
by 


c o)3 -1 o \ ^ i \ 


<1 ,2 1® 

^ -m6> j 

\3 


X _qf-.mca2y ^ Q 


(3.37) 
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^3.5 Discussion 

Experimental values of S are generally derived 
from the data on equilibrium concentration of defects as a 
function of temperature. The equilibrium concentration of 
defects of a given type p, at a temperature T is given by 



where /N, E is the energy required for formation of a defect 
and is the change in entropy (exclusive of configur- 

ational entropy) due to the defect. 

It has been observed experimentally that the specific 
heat of crystals near the melting points shows an anomalous 

increase from its expected linear behaviour due to the first- 

(3 2 ) 

order anharmonic terms . The anomaly in specific heat 

is attributed to the formation energies of vacancies. 

( 33 ) 

KraftmaKher et al. have reported a substantial amount 

of data on the concentration of vacancies in various 

substances (e.g, tungsten, molybdenum. Copper etc.) by 

measuring specific heat by modulation method (amplitude of 

temperature oscillations produced in the sample due to a 

current with constant and alternating components is related 

to the specific heat). Chekhovska and Zhukova have 

used the method of mixtures for concentration of vacancies 

in niobium. Calorimetric method has been used by Leibowitz 
( 35 ) 

et al. ' for measurement of. enthalpy of tungsten at 



37 


temperatures from 2800°K to 3600°K. Due to the volume of 
formation of vacancies, thermal expansion of solids also 
shows an anomalous increase at high temperatures. Kraftmakher^ ' 
has obtained a concentration of o. 8% of vacancies in platinum 
(near the melting point) from the anomaly in thermal expan'^^ 
sion. A discussion of such measurements has been given by 
Matthew (, The variations in A S of a vacancy^ obtained 
from experiments on different substances^ are rather large 
and it is inferred that some premelting phenomena and higher 
order enharmonic effects also contribute to the anomalies. 
MatthexAT has pointed out that these experiments might explain 
some premelting 'phenomena mainly as it is difficult to separate 
the effect of vacancy formation from other effects. 

A large number of accurate estimates of A S of a 
vacancy based on the lattice parameter and bulk density 
measurements have been reported. The difference between 
bulk density, p , and the density deduced from the lattice 
parameter, f^-^, is assumed to be due to vacancies such that 

% = (Px-P)/fx = , (3.39) 

where n is 'the number of vacancies in a crystal containing 

N lattice sites. A a and AL are the changes in average 

interatomic spacing a (obtained by X-ray method) and length 

L of the specimen, respectively. Recent papers by Smith 

and Chapman^ and Simmon et al. have listed a number 

(40) 

of experiments exploiting this method. Coufai et al. 
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have measured L of solid Krypton interferrometerically 

(accuracy more than 1%) for temperatures from 4 to 115°K. 

Interf errometric measurements have also been used for 

determining the equilibrium vacancy concentration in KCl 

(41) 

by von Guerard et al. 

For phenomena like diffusion, ionic conductivity etc. 
which depend on motion of defects, an equation similar to 
eqn. (3.38) holds for rate change. Measurements of ionic 
conductivity as a function of temperature gives inform- 
ation on /\S for a number of defects (log iT T vs. 1/T 
consists of straight lines in the high and low temperature 
regions due to different charge carriers). Experimental 

values are derived by least square computer fitting of the 

(42) 

data on . Allnatt and Pantelis have given a tentative 

explanation of the high temperature (>750° C) anomaly in ^ 

of NaCl in terms of cationic Frenkel disorder. Bauer and 

Whitmore have observed certain non-random deviations 

between the experimental results and theoretical results 

(based on presence of interacting cation and anion vacancies) 

on the ionic conductivity of NaP doped with CaF 2 which may 

be due to the presence of some other defects (e.g. tri vacancies# 

(44) 

interstitials etc.), Barr and Lidiard have reviewed the 
properties of defects in ionic lattices. 

Peterson^"^^^ has discussed diffusion in metals in detail 
in a recent review article. Mention may also be made of a 
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recent experiment on self -diffusion in selenium by Bratter 

and Gobrech^'^®^ which shows that the equibrium vacancy 

concentration in semiconductor's^ is several orders of 

magnitude smaller than that in metals. 

W, de Sorbo^^"^^ has give'n data on the energy evolved 

on annealing quenched in impurities in Au. Concentration 

of vacancies can be obtained by combining it with data on 

the change in resistivity due to formation of vacancies. 

( 48 ) 

Recently^ Polakova et al. have reported concentration of 

vacancies in nickel from quenching experiments. 

Among the theore-tical results^ the exact numerical value 
of AS of a substitutional isotope defect obtained by Nardelli 
and Tettamanzi ^ is based on the distribution function of 
square frequencies computed numerically by Overton and Dent. 
They have shown that the Einstein approximation and perturb- 
ation method give results consistent with their calculations 

only for very light isotopes and overestimate S o-therwise. 

( 14 ) 

Nardelli and Terzi find that the entropy of a vacancy 
at constant volume exceeds by nearly 10 percent the value 
obtained from the Einstein model. This shows, as pointed 
out by Nardelli and Terzi, ■t;hat Huntington et al.^^^^ have 
underestimated the contribution ^ to A S^of perturbation due 
to the loss of coupling around -the vacancy. Matthew^ has 
shown -that •the Einstein approximation gives exact results 
for isotope defects in cubic lattices ,.and^ in situations in 
which the force constants are changed^ it alwa 3 fs lies between 
■the exiact result and the trace approximation. 
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The pseudomolecular method for determination of 
As can be used for a general type of defect (i.e.,N’ 
may not be equal to N) by a proper choice of the pseudo- 
molecule, In the present model calculations of AS, the 
changes in force-constants due to relaxation of the medi\im 
around the defect can be included in the parameters occur- 

A. /V T A— 1 A 

ring in the effective force constant matrix^. A 

used for determining the pseudomolecular frequencies. 
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CHAPTER IV 


DETAILED DYNAMICS OF A POINT DEFECT 
^^4.1 Time-Dependent. Dynamics 

In this section we shall briefly review the use of 
partitioning technique for studying the time-dependence of 
atomic displacements in the defect region with certain 
initial conditions on the lattice. 

The equations of motion of a perfect lattice having 
N degress of freedom is (eqn. (2.1)), 


A 

M 


~ X(t) = L X (t) 
dt2 “ “ 


(4.1) 


From the theory of Fourier transforms it is known 

that if a function f and its derivatives Df, D2f, 

TC 1 £ • C K—1 

D ~-^f have jumps ^ o, ^o, ..,,o at t = 0 only, the 

Fourier-transf orm of the derivative of f, §[d^ f] , is 

given , 

J(D^) = (ito)^ §(f) _ |^(iw)K~l€Q + (ico)K-2 

(4.2) 

Taking the Fourier-transf orm of eqn, (4,1) and using 
eqn. (4.2) we get 


- u(co ) - j^i<M>X(o) -t X (O)] = li U(^ ) (4.3) 

where X(o) and X(o) represent the initial displacement 
and velocity vectors/ respectively. U(cp ) , the Fourier- 
transf orm of X(t), is defined as/ 
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U(co ) = ( X(t) dt 


(4.4) 


It is presumed that X and all its time-derivatives vanish 
for t < o. 

We partition the matrices of eqn, (4.3) as follows / 

/ uA^) ' 


'•A 


/ 

. / M 

M = i . I 

\ ° ™/ 




U(fu) = 


1 


(4.5) 


(L-M = 


A A O 


B- 


B 


; (L-Kca''^) 


_1 j Gg(0^ 

\ G0(<-3) 


The diagonal submatrices in eqns. (4,5) refer to 
regions 1*^2^ respectively, of the crystal. The 

region 1*^2 subsequently affected by the point defect. 

The matrices . corresponding to those defined in eqns, (4,5) 
for the imperfect crystal shall be primed except the sub- 
matrices and (Mj^) which remain unperturbed by the 

defect. 

Using eqns. (4.5) in eqn. (4,3) we get the following sets of 
equations , 


(A-f^^c->2)Uj^(tb ) + BU2(<^ ) = iC^]^AS3_(o) 4- Mj^^(o) (4.6) 

) + (C-m co^)U 2 (^ ) = i<i>iTK 2 (o) +n^ 2 ^G) (4,7) 

The required ec^ation for X 2 (t) can be obtained by 
eliminating U 3 _(<u ) from eqns, (4.6) and (4.7) and taking 
the inverse Fourier— transfoirm of U2^^ ^ obtained. The 
equation for U 2 (c->) is. 
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(C-W-B^(A-M^«2)“1S)u 2(«.) = i,jrfK2(0)-k-ffK2(0) 

- -1 (O) 

jr M^X^(o) J (4.8) 

We choose the initial conditions such that only the 
pseudomolecular atoms have non-zero initial displacements 
and all the atoms have zero initial velocity i.e.. 


2^1 <o) = 2i(o) = 2^2^°^ ~ ° 

and 

22(0) 0 


(4.9) 


This particular set of initial conditions is chosen as it 
simplifies the mathematics considerably and is applicable 
in many physical situations. 

Substituting eqns. (4.9) in eqn. (4,8) we get. 


Ox 'V / c- 2 . 

(C-m«"^)-B^(A-.M^o ) B 


U (<») = imWX (O) 
_2 —2 


( 4 . 10 ) 


Using eqn. (2,9) in eqn. (4.1o) we get the following 
equation for ^ 2 ^ ^ ^ 


U^C^) = i«:^G^(« )m X^CO) 


(4.11) 


In the presence of a point defect, equation for the 
Fourier transform of the displacement vector of the pseudo- 
molecular atoms, ^ ,can be written in analogy with 


equation (4,lo) as 



!ci 
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I C' ^B' I =iiTi' 


; ”2 


Using egns. (2.9) and (2.15) in eqn. 


r; (6J) = i « c ' -m’ <-0.^+ i («o) 


-C/ 


(4.12) we get 

g m*X' (o) 
aj -2 


(4.12) 


(4.13) 


4.2 Spectral Density Function of a Normal Mode 

A transformation to the normal mode co-ordinates can 
be performed by applying a unitary matrix S, constructed out 
of the symmetry adapted vectors of the point symmetry group 
of the crystal, to the displacement vector 
simplifies the calculations considerably as the various 
matrices occurring in eqns. (4.12) and (4.13) can be block- 
diagonalised by this transformation. 

Taking the Pour ier-trans form of eqn, (4.13) we get 


X' (t) 
“2 



j i4> I (O' -rh'ci^) +(G^^(<^) -64TTie^^) 
-oo 



m'X2^0) 


(4.14) 


Transformation to the normal co-ordinates can be written 
as 

Y^(t) = ‘§+X2^t) (4.15) 

Substituting eqn. (4.14) in eqn. (4.15), the -tilne-dependence 
of normal modes of the atoms constituting the imperfect 
pseudomolecule, Y^(t), can be written in matrix notation 
as 


Y' Ct) 
“2 


j J(«^)l2^0) 


io t. 
e d« 


(4.16) 
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where 

A 


5(«) = ^ \%+ g+^j. 


and 


Y'(0) = S'^X*(0) 
*“2 “"2 


(4.17) 


(4.18) 


From eqn. (4,16) the time-dependence of a particular 
normal co-ordinate (t) is expressible as/ 


yI it) = I J^('^ )e^‘''' "^dco , (4.19) 

— :0 

v/here 

J,, (^^) = I" iJi^)) (Y‘(0)) (4.20) 

^ ^ vyL. -2 

We use the following theorem to find the time-asymptotic 
form of y' (t) . It states that if a function f(x) satisfies 

i* 

the following conditions, 

(1) It has a finite number of singularities at X=X^, 

X2, X3, X^. . 

(2) f(X)-F^(X) has absolutely integrable N derivative 

in the interval including Xj^(for all m from 1 to M) , where 
Pj^(X) is a linear combination of functions having the 
nature of singularity at y;»Xj^,. ; 

(3) f^(X) , the derivative of f(X)^is well behaved at 
infinity, 

the Pourier-transform of f(x) viz. g(y) satisfies the 
following equation^ 
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g(Y) = + 0( {yT^) as |y| — / (4.21> 

m=l 

where is the Pourier-transform of P (X) . G (Y) for 

mm 

some of the functions Pj^(x)^ corresponding to some general 
types of singularities^ have been tabulated^ 

Using this result Montroll and Mazur and Baroody^^^^ 
have obtained the asymptotic forms of momentum autocorrelation 
functions which are written in terms of the frequency distri- 
bution function, Mon troll and Mazur have given results for 
the two-dimensional and three-dimensional lattices by finding 
the contributions from van Hove singularities of the frequency 
distribution functions. 

The exact evaluation of integrals in eqn, (4,19) is 
quite difficult because of the complicated nature of 
in most of the cases, i^) may be called the spectral 

density function of normal mode of vibration. We can 

study the asymptotic time-behaviour of yj_^ (t) by applying 
the above-mentioned theorem. Prom eqns. (4.19) and (4,21) 
it is clear that the isolated singularities and critical 
points of J-y)(co) will give non-zero contributions to 

3 ^ (t) for t Qualitatively, therefore, the well- 

behaved part of the spectral density function (cu) gives 
relaxing modes of vibrations and, the singularities and 
critical points correspond to persistent vibrations. Ihe 
time— dependence of a normal vibrational mode of the pseudo- 
molecule is thus expressible as a ;§uperposition of the 
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relaxing and persistent vibrations. The persistent modes 
having frequencies within the continuum shall represent 
resonance modes. An estimate of the width of any such 
mode, if present, can be obtained by substituting some 
values of the defect parameters and Green functions in 
the equation which determines frequency of the resonant 
mode. The singularities of J^, (co) outside the continuum 
give rise to localized vibrational modes. 

A sum-rule relating the trace of spectral density 

A 

matrix J(o>) of the perfect pseudomolecule with the 
frequency distribution function can be derived as follows: 
From eqn. (4,18) ) for the perfect pseudomolecule can 

be written as 

j(i.o) = iii s'* G (co)s m 
2 TT 

Trace J('^) = Trace G 

2Tr • c 

= - n giooo; t (4.22) 

2TT 

where the Green function g(000;‘^) is defined by 


g(000; ) 


i 2: i 

k,j 


(€ — ^o) 

(4.23) 


and n is the nurtiber of degrees of freedom associated with 
subspace of the perfect crystal. 

Prom equations (4,22) and (4.23) we have 
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Reel part j Trace J((U?) 


2"“ kTl 2 


Cl? n 
2 m N 




(4.24) 


2 

where p(^‘^ ) is the distribution function of squared- 
frequencies. 

4,3 Isotope Defect in the Linear Chain 

The time— dependence of displacement of an isotope defect 
of mass m' from its equilibrium position^ in a linear chain of 
atoms of mass m^and nearest neighbour force ‘.constant is given 
by using eqn. (4.14) as. 


'X(o) ^ g(0;^ ) 


x(t) = J f. , 2 77:-7V . (4.25) 

2 'ST 1 - ( m -m) C.O 2g(07 ) 

where X(o) is initial displacement of the impurity from its 
equilibrium position. 

Prom eqn. (4.25) 

- 'O 


X ( t) =- 


m'X(o) 

T\ 


o> Sin( <^t)g(0;^'^ 
l-(m' -m) &?'2g(o;a? ) 


_ 1 


Using -tiie principal part of g(o?<^) = ' ■ | a-^-, T o 

2rim o 1-Cos k)-«4 

in eqn. (4.26) we get 


(4.26) 




m 


iX> 


>^(0) Sin cot dc^ __ 

X V t) t X / o 2 ^ 


(4.27) 


where 


m 


>> = “Has and = II 

' m L m 
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For o < A < 


can be expanded as 


___ 

1 2 




(4.28) 


Substituting (4.28) in (4.27) we get 


X(t) = -nV ?- (l-A)"^ I (t) 

TT n=o n ^ 


(4.29) 


where 

I (t)= -A- 




<:o 


with 


"^n+l ‘^n,even') 

/ 


■^2 


S^_^^(t) = f Sin (n+l)9 Sin(63^ Sin© ) Cos 0 d(5 , 


'^/2 

C^_j,l(t)= J Cos (n+l)0 Cos( co^ Sin6 ) Cos© dP , 
and 

j (co t) is a Bessel function of the first kind, 
n h 

This example clearly shows the difficulties involved 
in determining X(t) explicitely for a general defect. 

Takeno and Hori^^^^ have shown that a very heavy impurity 
in one-dimensional and three-dimensional lattices behaves like 
a Brownian particle in the continuum approximation. Recently, 
Eisenriegler^S'^^' has discussed the dynamics of an isotope 

2 r 2 ~2* • 

defect analytically assuming a' spectrum of^o^ -<p 

the frequencies of the perfect lattice. 
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4 . P'^'naraiGg of the Atoms around a Vacancy 

Vacancies in a lattice have received particular attention 
in literature because of the role they play in raany physical 
Phenomena and the vast amount of data available on them. However^ 
no exact method of studying the dynamics of a vacancy seems to 
have been given so far. The model of a vacancy that is generally 
used is to regard it as a limiting case of a substitutional 
defect of mass m' and force constant qp' for m' 0 and 
1’ 0- Effectively, this defect is know to weaken the 

force constants considerably. In an alternative model, therefore, 
vacancy at a lattice site is produced by cutting off the force 
constants connecting the atom at that site to other atoms; mass 
of the atoin at the vacant site being immaterial. 

We shall briefly mention seme of the results about vacancies 

( 18 ) 

obtained by earlier authors. Lengler and Ludwig have shown 
that a vacancy in the simple cubic lattice does not have localized 
modes of vibration. This result has also been reported by Janot 
and Champier^^^\ 

Land and Goodman^ have used a molecular model to 
investigate the possibilities of localized modes at a vacancy. 

They assume, to begin with, that only one or two shells of atoms 
around the vacancy can relax and calculate the frequencies of the 
nommal modes. The highest! freqiency modes obtained from them are 
then identified as the possible localized modes,, the changed 
frequencies of such modes, ‘with increasing number of free shells 
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around the vacancy are then calculated. Land and Goodman 
iicivc proven tsiat the existence of localized modes at a vacant 
si .e deoenas r>.i3.inly on the relaxation of structure when the 
atom is removed and consecpjLently on the form of interatomic 
po"'ential and lattice tyt>e. A vacancy, for examole, has a 
localised mode (v/ith frequency 1 % above the continuum) of ^ 2 ^ 
syiccetry in copper but does not have localized modes in iron 


and sodium, Takeno 


( 57 ) 


has discussed the vibrations of a diatomic 


simole cubic lattice v/ith a hole. He has ignored the changes in 

■oola’^izabilitics and effective charges of atoms and the vacancy 

has been regarded as a change in force constants only, in the firs- 

approximation. Takeno has shown that the vacancy may have 

localized modes in the forbidden gap but has no localized modes 

outside the optical band. 

( 22 ) 

Krumhansl ' has shown that lattice vacancies can act as 

resonant scatterers at frequencies well below the Debye 

frequencies. Such a resonance scattering has been observed 

experimentally by Schv/artz and Walker They have measured 

the thermal conductivity of alkali halides doped with divalent 

impurities. Schwartz and Walker have interpreted the dip in 

the thermal conductivity curve as mainly due to the phonon 

resonance as it is found that the position of dip is independent 

of the mass of the divalent ion and^ its agreement with the 

( 22 ) 

theoretical results of Frumhansl' requires practically no 
force constant for the scattering centres. 

' LI.T. KANPUR, ' 

I —J TWAU LWW AHT, i . 
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4 , 4a Vibrational Mod^s of a Vacancy In the Square Lattice- 
A vacancy is created at the origin of subspace 
cutting off the force constants connecting the atom O to its 
four nearest neighbours labelled as 1 to 4 (figure (la)). 

The atoms (1—4) constitute the pseudomolecule. In the scalar 
model/ the displacement vector X^^t) of the pseudomolecule is# 


X' (t) 
“2 


XjCt) 
XjCt) 
X4(t) 


(4.30) 


where is displacement of the n^ atom from its equil- 

ibrium position. 

The unitary matrix s' representing a transformation to the 
no 3 rmal co-ordinates of the pseudomolecule is given as 


S 


1/2 

1/2 

l/f2 

0 \ 

1/2 

-1/2 

0 


1/2 

1/2 

-1/72 

0 

1/2 

-1/2 

0 



(4.31) 


G^(co) and 8(^3, 4a)/ it can be shown 

fD-4'?r 4 tooj2 E 


G ) — C— in 

c 


E 

P 

E 


D-4'T 

E 

F 


F 

E 

D-4T +m 
E 


E 

P 

E 

D,4T -hncc' 


(4,32) 
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v;here 


^ + 5<9o -^2-2311) + (3o-32>"^ 


1 T Oo 1 , . ^ .“1 

5 5 51 ■ 2 + 92-29ii> 


"""* “ 2 R I " 2 ‘®o^2-29i1>‘' - 


(4.33) 


The spectral density matrix of normal modes of the vacancy 


J(<^) = 


/-(Ho 

gi 


(go+g2"^^ii^ ° 


(go-g2^ 


(go-ga^” 


' (4.34) 

Prom the discussion of ^4,2 and eqn, (4.34) it follows that 
the frequencies of persistent modes^ if present, are given 
by solutions of the following equations (the equations for 

o 2 

/ Co arise by equating the real parts of basic equations 
^ L 

to zero to look for possible resonances) . 


(i) Mode (1,1, 1,1) : go = g^ (4.35) 

( 3) 

Substituting the analytic forms of the Green functions' 
and the sum-rules satisfied by Green fimctions (eqns.'(3.9) 
and (3.1o)), eqn. (4.35) gives the following conditions: 

KCVj^) = - , 0<‘02<;l.a2 (4.36a) 
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1 + >1 



< Co2 < 6 o2 


(4.36b) 


K( p ) 


where 


■n = 

P = 


4?'-nT^ 

4 y 

__4^__ 


■' ( 1 ) 

2 4+f ' 

2 


^ I 

-1 


2 2 
CO 


mco^_4 T 


b m 


4Tr 

8T 


(4.36g) 


and K's are the complete elliptic integrals of the first 
hind. None of the equations (4. 36a-4. 36c) has a solution 
(Figure 3a) . 

(ii) Mode (1, 1,1 J) s ^ 

Eqn. (4,37) yields 

( 1 - /|^)K()^ )=2E( Tf|)- ^(1+1^), 0 <c0 < |ch2 (4.37a) 

^ 2 

(1- 7| * ^) K( V|^ '.) =2E( ■Yj ‘ ) + -—(1- 'T| ' ) , ^ (4.37b) 

(1~ p^)K(p )=2E( p)- ^(1-f ), ci (4.37c) 

where E ’ s are the complete elliptic integrals of second kind, 

Eqns. (4, 37a-4. 37c) do not have any solutions (Figure 3b) 
showing thereby that (1/1/ 1,1) is not a persistent mode. 

(iii) The doubly-degenerate mode (l/O/l/O) and (0, 1,0, 1) : gQ~g2“^4^ 

(4.38) 


Eqn. (4,38) leads to 


2 

( 1 _ 'y|^)K('>| ) =E(o^ )-► ^‘(1-11 ) , 0 < cu^ <. (4,38a) 

(l-ij ’^) kC '»|')-E( 'VI’)- ^{1+1|’), (4.38b) 


eC p) 


it(l+f ) 


2 ^ 

to < ^ 


(4.38c) 
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Eqns. (4. 38a-4. 38c) also have no solutions ^Figures 3a & 3b) 
Thus we conclude that a vacancy in the square lattice can 
have only relaxing normal vibrational modes. 


^ 4 . 4b Vacancy in the glmple Cubic Lattice 

The model used here is the same as that discussed 
in Chapter III. The displacement vector of the pseudo- 
molecular atoms is (Figure / 

/xi(t) \ 


x'Ct) 


Xn (t) 


Xj ( t) 


X^(t) 


XgCt) 


Xg(t) 


(4.39) 


vrhere Xy^('t) is the displacement of atom. 

The matrix representing transformation to the normal coordinat 


Is given by 

i 

1 i/F 

1/J2 

0 

0 

l/JlT 

1/2 

1 

1/16 

0 

1/.J2 

0 

i/JiT 

-1/2 


1/J6 

0 

0 

l/^'2 

- 2 /JT 2 

0 


1/J6 

-i/i'z 

0 

0 

l/.p^ 

1/2 


1 1/^ 

0 

-1/.J2 

o 

1 /J 12 

-1/2 


\iAl^ 

0 

0 

-1//2 

- 2 /JI 2 

0 





(4. 

40) 

The labels of the 

rows and columns 

of in 

eqn. (4,4o) 



are in the order of indices of atoms 1-6, 

The form of spectral density matrix J(co ) is given 
by 
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/ g 

-1)^ 0 


) = 




2T?i 


0 

0 

0 


0 


(g„-g2) 


-1 


- T 

f ' 

0 


u 


(go-g2^ 


-1 


0 

o 

o 


- Y 

0 




-1 


0 

0 

0 

O 


- T 
0 


(go+g2~ _i 

-2g3) 

-T 

0 


o 

0 

0 

0 




-2c 


with cu) = g(4-.&[, V-^'^--C3-5^7CD) . 


(4.41) 


= g( 0 , 0 / 0 ; co) , = g(l, 0 , 0 ;co) ; g 2 =g( 2 , 0 , 0 ; c^) 

and g 3 =g(l/ 1^0; <g) and these elements of the Green functions 
are as defined by eqn. (3.21) ■ 

The conditions for existence of the persistent modes are/ 


(1) 

Mode (1/ 1/ 1/ 1, 1/ 1} 

o^ 

1! 

(4.42) 

(2) 

Mode (l,O/0/l’/0/0) 


(4.43) 

(3) 

Doubly-degenerate mode 

(O/l/O/O/l/O) / (O/O/ 

l,0/0/l) : 


To find the solution of these equations the required 

Green functions have to be computed numerically. The method 

of Fourier-Series expansion (viz, the imaginary part of a 

Green function is expressible as a Fourier sine-series with 

some Bessel functions as coefficients) has been used for the 

same. It is found that only the mode (1/ 1,-!, 1/ 1, 1) has a 

2 2 

solution which is a rather broad resonance at co = o.S^^L 
(Pig. 4), The other two equations, (4,43) and (4.44) do not 



Re g^Co)) 



Resonance mode of a vacancy in the simple cubic lattice-. Root of Eqn Re 90(00) = 
Re 04(00), 
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have any solutions and therefore corresoond to relaxing 
modes of vibration. Thus a vacancy in the simple cubic 
lattice does not have localized modes of vibr-^tion. 

The result that a vacancy does not have localized 
vibrational modes is understandable on the basis of a 
substantial amount of vreakening of the force constants 
caused by it. By Rayleigh’s theorem^ the frequencies 
of such a defect should decrease. However, relaxation 
around the vacant site may adjust the force constants such 
that it shows localized modes, as has been discussed by Land 
and Goodman Many authors have given methods of finding 
relaxation around the defects in cubic lattices. In more 
detailed calculations , relaxation due to defects can be 
taken into account by choosing a bigger oseudomolecule. 

4 . 5 Correlation Functions Associated with a Point Defect 
The formalism of ^4.1 and ^4.2 shall be used in this 
section to give a method of comouting the correlation 
functions associated with a general ooint defect which are 
required for calculations of Mdssbauer line-width, neutron 
scattering cross-section, optical absorption, transport 
Coefficients etc. . 

The time-dependent correlation functions of displac- 
ements and velocities of atoms in the imperfect crystal 
( 59 ) 


are defined as 



6t 


4- os 

<x (t7t)x (t‘j(0)> = i ( dcoexp(-icot) f g - — 

' ^ f ^ exp(-p^)-l 


Im G • ( 4, 1‘ ; co^-ie ) e— ; o 




(4.45) 


> t) p U * ; 0)^ = 


ti 

— MM, 

rr 4 i 


.£30 


( dco.exp(-iost) 

exp(-p]lo)-l 


Im G'^ U,V; c3-i&) 


(4.46) 


^X^ (£_; t) p^ (V ; 0^ 


m 


-r 

M^, J d<s-> exp(-iwt) 


Co Sgn 


■g J ■ expC-pti^^) ~1 

Irn (i,J'7o-ifc), (4.47) 

where /(■ ■ - "^denotes the thermal average; = l/3cgT, kg 
being the Boltzmann’s constant and T is the temperature, 
b*, the Green functions matrix of the imperfect cr^^-stal, is 
defined as 

b' = (L'-M'co^)”^ (4.48) 


where L’ and are matrices of the atomic force constants 
and masses , respectively. 

For correlation functions associated with a point 
defect^^and i.’ in eqns. (4.45) -(4. 47) are to be restricted 
over the indices of pseudomolecular coordinates. Therefore 

A ^ 

G' in these equations can be replaced by the stibmatrix G^ which 
can be written in terms of G^ and parameters of the defect 
by using eqns. (2.9) and (2,15). Hence eqns, (4,45)-(4,47) 
can be revtitten as. 




ft 

TT 


/ dcj, 


•iot Sgn 




exp(4M-i 
—1 .-i A 2^ "7 


i 0^ 


c,a 


Im 


(4.45a) 
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JfCO 


< P/i/- 1) p (Si ’ ; o)> = 2L M M . r do e-i ■t -.Ci£sgn.,id. 
^ f>" ^ TT -e -e -L^ e - _ 


Im 


0 ^ — 1 A 9 


-1 -1 


(4.46a) 


4.4?o 


<x(£;t)p(l';0)> = is M , r ^-.■- 


A . 2 9 ^. 


rn'c.7^-C' +(G^‘-lC-mao2) 1 

c ^ a J 


Im I m'c. 7 ^-C' +(G^ -iC-fnao^) | (4.47a) 

Mansel et have studied the Debye-Waller factor of 

the Mossbauer isotope obtained by irradiating platinum wi 

neutrons. They have concluded that the Debye Waller factor dro 
appreciably when the Mossbauer atom is nearest to a close Frenk; 
pair. As an example, therefore, we shall obtain here an express 
for the Debye Waller factor, f, of one of the nearest neighbour! 
of a vacancy in the isotropic simple cubic lattice (model consi< 


ed here is the same as that used in ^4. 4b). f is defined as^^^, 


f = exp (-2K<(_ x2( lo)^ ) 


(4.48) 


<^x‘^(lO)^is obtained from eqn. (4.45a) by putting t = 0 


.CO 


<^x^(lO)/ = j dco Goth ( Im G'(lo,10; cjo -ife) (4,49) 

Using the matrices required in eqn. (4.49) from those 
given in ^4, 4b we get 


Im G^(10;10; 0^-14) = i Im 


Pq. 


2(go+g2"^^3^ 


(Po-gi) 1- ■s'(go+92"2g2^ 


3 ( 90 - 92 ) 


l- TT ( 90-92) 


] 


(4.50 


Imaginary parts of the various terms occurring in eqn 
(4.50) can be simplified as. 
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^m{g^-g2) 

( 1- YRe(go-g2) ) 2+(im(go-g2) ) ^ 


( 4 . 51 ) 

Ihe required Greens functions can be computed numerically . 

We have not cojrfcuted any numerical results here. 

Lakatos and Kruirihansl^^^^ have given analytic expressions 
for the mean square displacement of the impurity and its 
nearest neighbours in f.c.c. and b.c.c. lattices. Goldman^ ' 
and, Haridasan and Nandini^^^^ have reported calculation of 


Debye-Waller factor in rare gas solids and cesium-halides. 

Calculations and precision measurements of Mossbauer fraction 

of Fe^”^ impurity in Pt, Pd and Cu have been reported, 

. (6 

respectively, by Patnaik and Mahanty^®'^' and Nussbaum et al. 


64 


CHAPTER V 

PSEUDOMOLECULAR MODEL FOR A SUBSTITUTIONAL IMPURITY 
IN A LATTICE OP POLARISABLE ATOMS 

^5,1 Introduction 

In the calculations based on pseudomolecular model 
in the previous chapters, the atoms were regarded as non- 
polarizable. The usefulness of the shell-model treatment 
of lattice vibrations has been amply demonstrated by 
now^^'^"^ , In this model, as proposed originally by Dick 
and Overhauser , each atom is regarded as consisting of 
a core surrounded by a massless shell composed of some of 
the outermost electrons. The cores and shells are assumed 
to be connected by certain short-range force constants. 
Under the influence of these short-range forces, each atom 
has ■ a distortion polaris ability in addition to the polar- 
izability induced by the electric fields. The potential 
energy for such a model, in the harmonic approximation, is 
expressible as a quadratic in the core-displacements and 
electronic dipole moments. In this chapter we shall cons- 
ider a pseudomolecular model for a substitutional defect 
which does not change the long-range Coulomb forces due 
to the charges of shells and cores. 

{ 69 ) 

In the notation of Cowley et al. , the equations 
of motion of the lattice in the she 11 -model can be 


written as. 
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/ 

V 


oc - «c"' 


6'" 

-cs 


^cs 


4; 


ss / 


u 


I ^ / 


= 0 


(5.1) 


where U and V denote the displacements of the cores and 
shells, respectively, from their equilibrium positions, 
^cc' ^v'’" ^ matrices of the force constants 


:cs ^ss 

betx-^een cores and cores, cores and shells and between 

shells and shells, respectively. Each of these <p matri- 

R 

ces consists of a short range part and a long-range 


part (arising from the Coulomb forces) such that 

'5’gc 


tR A ^ 

_ + X C X ; 


and 


I 

^SS 


'•cc 


^ ^ A A 

(±> + Y C Y , 
^SS ' 


I 


cs 


(^^) 


/\ R A A 

= . + X C Y 


cs 


(5.2) 


where X and Y are diagonal matrices of the core and shell 

A 

charges, respectively, C is a matrix of the Coulorrib coef- 
ficients. 

Eliminating V from equations (1) we get 


—1 

' - icsiss %s'‘~ = ° 


( 5 . 3 ) 


-1 

^R ^ R 
^ss 4 cs 


is a long-range 


(70) qiR 

As discussed by Lax , 

term due to the occurrence of ('$ss' Though, in 

principle^ eqn. (3) can be partitioned in the same way as 

done in Chapter II, this is of no special advantage here as 

the long-range nature of %g 4 ss “^cs cumbersome 

to take into account changes in the shell parameters at the 

f 
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point defect. We shall show here that if the core and 
shell coordinates are taken together to form the pseudo- 
molecule and if the long-range Coulomb terms in eqns. (5,2) 
and (5.3) are assumed to remain unchanged, a local effective 
interaction term can be obtained in terms of the defect 
parameters . 


5,2 Pseudomolecular Model 


We partition 


U 

V 


and the matrix of force constants 


in eqn. (1) as follows 


U \ 


\^i 


21 

22 

I 


(5.4) 



n. 


^2 


n. 



T- 9*^ X 

^ $cc“ ^ 11 


'4s>ix 

^ ^cs^ 12 1 

^ lcc“ 

( $cc“ ^ 22 

^ $cs^ 21 

! 

^ 'fcs' 22 

^-Ics^ll 

^ ^cs^ 12 

^ 4s^ 11 

'4s>12' 




/V 

^ ^cs^ 21 

^ 22 

^ iss^ 21 



(5.5) 


where 1 and 2 refer to indices of regions I and II of the 
crystal, n^ and n 2 are degrees of freedom associated with 
regions I and II, respectively. The dimensions .of-v 2 . --f 
the submatrices of ^ in eqn. (5.5) are as indicated by n*s. 
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We define a unitary transformation matrix S as 
follows ! 


s = 


ni 


Hi 


n. 


n. 


n. 


I <i- 

1 I 


0 


^^2 

0 

0 


0 


n. 


0 


/V 

I 


0 


n. 


O \ 


O 


(5.6) 


Applying the transformation § bo eqn. (5,1), with the 
displacement vector and matrix of force constants parti- 
tioned as shown by eqns. (5,4) and (5.5), we get 



(5.7) 

(5.8) 

(5.9) 


From eqns. (5.7) and (5.9), eqn. for can be written 
as 

or 

X 2 ° • (5.11) 
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where 


I $ 


s \ 
ff- \ 


-S i '4 , 


11 't 12 


t .h® 


"^12 


■^22/ 


From eqn. (5.8) 


and G„ 


i ^ A 

?-i 1°^ 

^3 = e, 

\ 




r“l -T- ^-1 G-1 
*^3 “ 'is ~ ^ ^ ^ 


A, /S. 1 

S G S 


(5.12^ 


(5.13) 


and 


'''t - .s2m 

±CC 60*’^ 


r~l 

G = <^ JT 

\ ior. 
-\ 


h H -1 
^cs \ 

Issi 


‘ G G 
CC CS 


8”^ G 
CS ss 


(5.14) 


The required submatrix 8^, in eqn. (5.11) can be 
shown^ from eqns. (5.13) and (5.14), to be obtainable from 
the submatrices of G as follows: 


/ (Sco) 


cc^ 22 


A f 

■ = -r 

\ 22 


2Z \ 
22 / 


(5.15) 


where indices in the submatrices of G^ are restricted to 
those of region II only. 

The matrix of Green functions, defined by eqn. (5,14), 
can be written as 


2n- 

Scc(S-B‘;“) =-^ e 


-ik, (R’-R) 


D (k) 


(5.16) 
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Gcs(R,R';co) 



dk ®’-£) 


D 

CS 


(k) 


D3g(k)~Des(k) 6^3 (k) 


(5.17) 




2^ 


Ggg (R 


= -TT^ r 

2fT j 


dk e 


-ik. (R'-R) 


(k) 


ti^c(5^)6ss(k)-1ics(k>®cs^^> 


(5.18) 


where 


°ss‘S> ' $ss'£'S'> e"--- 


los'S.S') (5.19) 


and 




2^ X ik.(R‘-R). 


cc - R/R‘ 


We shall consider here only those substitutional 
impurities which do not change the nui^er of degrees of freedom 
associated with the lattice .and only the core and shell 
parameters of the atom replaced by the impurity are affected. 
The matrix C in eqn. (5.2) is also assumed to remain unchanged. 

The formalism so developed may be useful in many physical 

*1^ 

situations (e.g. # for defects like Li and K in Nal crystals 

(71) 

which have been treated as essentially mass defects j. For 
such defects, the submatrices and of cp® 

(eqn, 5.12) remain unaltered if the pseudomolecule (i,e, ^ 
is chosen such that the short-range shell parameters between 
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the peripheral atoms of the pseudomolecule and those of 
the unperturbed region I do not change. The change in 
matrix 0 ^^ *^t.e to the defect can be determined in terms 
of the shell-'^arameters of the perturbed region II only. 

The equation of motion of the imperfect pseudomolecule 
can therefore be written from eqn. (5.lo) as/ 



(5.20) 


where prime denotes the perturbed matrix. 


5,3 Application to Simple Cubic Lattice 

As an illustrative example we shall obtain here the 

A I /V 

form of the perturbation matrix <i>® for a substi- 

^22 “ 22 

tutional impurity in the isotropic simple cubic lattice. In 

A g 

the nearest neighbours interaction approximation, $ and 
^ s * 

^ are (14x14) matrices. Various shell parameters between 
2 2 

any two neighbouring atoms of the perfect lattice are as 
depicted in Fig. ( 5 ) . 



2 ( 010 ) 



Shell model parameters for the pseudomolecule of the perfect simple 
cubic lattice ( § 5-3). 
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Time-independent equations of motion of the core and 
shell of the atom can be written as, 

-Me Oih = d( 2Ui- 6U )+F( lv.-6U )+K(V -U_) , (5.21) 

i=nearest i=nearest 

neighbours neighbours 
of n of n 


0 = N ( IVi-6Vf^) +K(U^-Vj^) +P ( Ui-6V^) , (5.22) 

i=nearest i=nearest 

neighbours neighbours 

of n of n 

where is the mass of the core of the nth atom and shells 
are assumed to be massless. 

Matrices 4' cT; and defined by eqns. (5,1) 

and (5,2) can be obtained from eqns. (5.21) and (5,22) written 
in matrix form. Using them in eqns. (5.8) and (5.12), matrix 

e ^ S * 

cis^ can be obtained. Similarly, ^ matrix can be obtained 
^■22 -22 

by xvriting the equations of motion for the imperfect crystals. 

The explicit form of ^ terms of the changes AN, 

22 ^ 22 

.’iF, d and in the shell model parameters due to the 
impurity is given in Appendix VT . 


The procedure for obtaining the spectral density function 
associated with each mode is similar to that given in ^4. 3, 

The eigen-value equation obtained for the doubly-degenerate 
E,.-modes (0112 112 OOOOOOO) and (OOOOOOO 112 11*2) , for 
example, is 
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/ _ G (o)-2G (3)+G (2) G (o)-2G (3)+G„ (2) \ 

j Ad+ ^ -p.;. ss ss ss cs ^ cs \ 


^ 1 At <V 


^SC (0)-2Gs;,{3)-W3g^(2) 


^SC 

^ - I . _. 

'^Go'^Ie, 


'^Go'^Ie 


AN+ AP+ 


g 


*^cc (0)-2G^^(3)+G^p(2) 




cc 


CC ' 


V^G^V 




where 

A-l-A A 


Ul-Hi2’-2u3HfU4+U5-2u5 


t 

i=0 


|v^+V2-2v3+V4+V5-2vg| 


(5.23) 


= '°cc(O)-2Gcc<3>«oc(2>)tG^^(0)-2G„,(3)-K3„(2)) 


SS 


ss 


-«2sc 'O) -2G^^ ( 3) G3^ ( 2) ) (G^3 (0) -2G^„ ( 3) -K3^, ( 2) ) , 


CS 


cs 


(5.24) 


and 


= G( R-R‘ ) = G(R.-R', R-,-R'/ R^-R') 

i-j_223 3' 

G(3) = G(1,1,0) 

G(2) = G(2,0,0) (5.25) 

G(o) = G(0/0/0) 

G(l) = G(1,0,0) 

and these elements of the Green functions are defined by 
equation (5.16-5.19). 


5.4 Conclusion 

It may be mentioned here that the formalism in §5,2 

( 19) 

is quite similar to that given by Page and Strauch 
using Lifshitz's approach. Their method enables one to 
take into account the changes in Coulomb terms due to the 
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changes in core and shell charges (changes in the Coulomb 

matrix due to relaxation of the medium around impurity 

have been neglected) as a local non-zero submatrix in the 

full perturbation matrix. The dimension of the force 

constant matrix in their treatment is increased due to 

an additional (linearly deoendent) equation of motion in 

eqn. (5.1) and calculations of the Green functions are 

therefore more involved -than in the present approach. It 

( 72 ) 

has been pointed out recently by Schroder that a 
'Breathing shell model' , in which the electronic shells are 
regarded non-rigid, gives a better agreement with the 
experimental dispersion curves of some alkali halides. The 
present method and that suggested by Page and Strauch can 
both be generalised to the 'Breathing shell model*. 
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CHAPTER VI 

CONCLUSION 

The pseudomolec\alar model of a point defect used 

in this work is a development on the concepts used earlier 

(7 8 9 29 ) 

by some authors * • ' ^ 1^3 applicability to the self— 

entropy problem has been discussed by Mahanty^^^^. The 
advantage of this method over the earlier exact and appro- 
ximation methods lies in the fact that it is readily applic- 
able to a general point defect which may involve a change 
in the number of degrees of freedom associated with the 
lattice. Also it is an exact method with computational 
problems not more involved than what is required in conven- 
tional methods. This work contains only model calculations 
for defects in square, simple cxibic and b.c.c. lattices. 
Calculations should be done for point defects in f.c.c, 
lattices for v^ich a lot of data exists. 

Neglecting relaxation of the medium produced by 
point defects, the pseudomolecular model has been used to 
study the exact detailed dynamics of general point defects. 
The method leads to the exact result that a vacancy in the 
simple cubic lattice has only relaxing vibrational modes. 
Land and Goodman^ have pointed out that localized modes 
of vacancies (in cubic lattices) with frequencies slightly 
above the continuum are possible if relaxation of the stru- 
cture is taken into account. Relaxation of the medium can 
be taken into account by choosing a bigger pseudomolecule. 



A sum-rule for the spectral density matrix of the perfect 
pseudomolecule has been given. Some similar sum— rule for 
the imperfect pseudomolecule, if possible, may give inform- 
ation about the frequency distribution function of the 
imperfect crystal. 

The phonon-scattering from point defects has not been 
discussed within the framework of this formalism. Using 
a slightly different partitioning method the phonon-scat- 
tering can be described in terms of the usual t-matrix 
formalism, as has been discussed in detail by Wagner . 

The correlation functions associated with a general point 
defect, required for many physical problems e.g, scattering 
of particles by the lattice, Mossbauer effect, optical 
absorption etc., can also be obtained by the pseudomolecul- 
ar model. In the present work, even though results for 
practical models have not been computed, the Debye-Waller 
factor of one of the nearest neighbours of a vacancy in 
the simple cubic lattice has been given in terms of the 
parameters of the lattice. 

The theory of pseudomolecular model for the substitution- 
al point defects involving a change in the short-range forc- 
es in lattices with polarisable atoms has also been given. 

The possibility of application of such an approach to 
general point defects which may change the long-range 
Coulomb interactions has not been investigated in the 
present work* Such an extension may provide a very \iseful 
model for the colour centers in ionic crystals. 
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APPENDIX I 


The block— diagonalised form of (eqn. 3,3) for a 

substitutional defect in the simple cubic lattice 



with 


= § ( 1 -S) ^ 


^2 " -2 o:^G( 100,000) 


^^5 T 2 J 


r m 

^/AT2 "^2 

r 


A3 = ^(-4^) cu^G (001, 000) r ^ 1 

M 1 AT2 'T2 J 

= 2(14‘) G(000,000)-G(100,000)J 


B 2 - 1- 


M 


B3 = 2/2 ci~S) 


J^4G ( 100/000) -G ( 100 / 100) -2G ( lOO, OlO) -G (OOO, 000^ 
[g { OOl / 000 ) -G ( 010/ OOllj 


T 2 

M 
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= 2 j2<K(l~S) G(010,000)-G(010,001)J 

= 1 + 'X (1-5) •— .‘2G (001,000) -G (000/ OOO) -G(00l,00l^ 

/K'^o I — 

= 1 + -=^ |G(000,000)-G(100,100)| 

= 1 - [2G( 100,010) -G(l00,l00) -G(000/ 000) I 


1 +rA (000,000) -G( 00 1,001) 


where 


and ^2 the central and aon-central nearest 


neighbour force constants ( =c< ) 


7r2 +A32 


m = mass of a normal atom. 


m+ (5m and T + Al are the mass and force constant of the 

substitutional impurity. 

and 

A , E / B and A. are the irreducible represent- 

Ig t ig 2u 

ations of the point symmetry group 
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APPENDIX II 


G^{c^) for the Isotropic Square Lattice (^3. 4a) 


0 


0 jg( 0 ,Q 70 ^) g(l,0;<:^) g(l, 0 ;CJ) g(l,0;<^' 

g(l,0?^) g(0,0;t^) g(l,l;c.o) g( 2 , 0 ;co) g(l,l;cu) 


A 

Gc(<^'') = 2 I g(l,0;w) g(l,l;c^) g(o,0;uj) g(l,l;U) gi2,0;c^) 

3 lg(l/0;c.^) g(l,0;<j) g(l,l;4-^ g(0,0;f^’) g(l,l;cO) 

4 W(l/O;<o) g(l,l;<.-J) g(2,0;‘^) g(l,l,*<o) g(0/0;<*j) 


where 


(o, 0) =G^ ( 1 , 1 ) =G^ ( 2, 2) =0^ ( 4 , 4 ) =g(o, 0;M ; 

G^^ ( 1,0) =Gc(0, 1) =G^ (2,0) =G^(0, 2) =G^ ( 3,0) 

~Gc ( 0, 3 ) “G^ ( 4 , 0 ) ( 0 / 4 ) —g ( 1 , Ojco) • 

G^ ( 2, 1) =G^ ( 2) =G^ ( 2, 3) =G^ ( 3, 2) =G^ ( 3,4) 

=G^(4,3)=Gc(l,4)=G^(4,l)=g(l,l;co); 

and 


Gc(l,3)=G^(3,l)=G^(2,4)=Gc(4,2)= g(2,0;^) . 
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appendix III 

hatrices S and V for the isotropic square lattice (^3. 4a) 

4 

0 

0 


S 


A 

V 



0 

1 

2 

3 

0 

) 

, 1 

> 

0 

0 

0 

1 

: 0 

f 

i 

1/2 

1/2 

l/f2 

2 

0 

1 

! 

1/2 

-1/2 

0 

3 


1/2 

1/2 

-1/J2 

4 

\-° 

1/2 

-1/2 

0 


0 

1 

2 

3 

° 1 

1 Vfs 

-4/ 20 

0 

0 

1 

l/Js 

1/ 20 

1/2 

i/n 

2 

1/45 

1/ 20 

-1/2 

0 


l/l5 

1/ 20 

1/2 

- 1/42 

} 

4 

\ ^ 
U/s/5 

1/ 20 

-1/2 

0 


1/J2 


■ 1 /^, 


4 


0 


1/42 


-l/f2 



o> 
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APPENDIX IV 


r- 

V 


= 3 


the 

isotropic 

simple 

cubic 

lattice (^j3.4b) 

0 

1 

2 

3 

4 5 

6 

l/-i7 

6/U2 

0 

0 

0 0 

° \ 

1//7 

-1//42 

1/72 

0 

0 l//r2 

1/2 

1/77 

-1/742 

0 

1/72' 

0 1/712 

- 1/2 

l/JT 

-1/742 

0 

0 

1/72 - 2 //i 2 

0 

1/77 

-1/742 

-l/v/2 

0 

0 i/7i^ 

1/2 

1/77 

-1/7^ 

0 

-1/72 

0 1/712 - 

-1/2 

\l//7 

-1/7^ 

0 

0 

-1/72 - 2 /J 12 

0 . 


'6T -7 
-t 6T 

-If 0 

-t 0 

-Y 0 

-Y O 

-Y O 


-Y 


0 


6 Y 


0 


O 


O 


0 


0 0 

0 0 

6Y O 

O 6Y 

O 0 

O O 


-Y -T 

0 0 

0 0 

0 O 

O 0 

6Y 0 

O 6T 
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/^+f. 


j 

-5g^) 


+4g-3+g ) 
z 


3 

-5g^) 


0 

0 

0 

0 


i(7g^+g2 


7 


y3-l2g^) 


+4 a — 
0 
0 
0 
0 


0 


(Oo-gj) 

o 

o 

o 


0 


0 


0 

(go-g2^ 

0 

0 


0 

o 




0 


o 

0 

o 

0 

o 

(g„+g,) 

'■u Z 

-Zgj) 


'3o- 

-2g. 


where 


3 A-1b = 


/ 0 
0 
0 
0 
0 
0 




o 

p 

Q 

Q 

R 

Q 

Q 


0 

Q 

P 

Q 

Q 

R 

Q 


0 

Q 

Q 

P 

Q 

Q 

R 


C 

R 

Q 

Q 

P 

Q 

Q 


P = 6T-[| + Ug^-g^r^ + l(g^+gj-2g3)‘^^ 


° = I ■ I ‘9o'*^2-'33> 


0 

Q 

R 

Q 

Q 

P 

Q 


and 


R =t 


1 (g„-g,)“^ + i (g^+gj-sg,) 
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with 


/V. 

V 


11 


appendix V 

V I-IATRIX FOR THE B.C.C. LATTICE 


/ 

/ V 


V 


11 


A 

V 


V 

A 

V 


21 


31 


-8/J72 O 

i/fn 0 

lAi72 O 

i/fn o 

1//72 0 

1//7? O 

i/fn o 
l/fi2 0 
1/m o 


V 


V, 


o 

0 

0 

0 


12 

22 

32 

0 

0 

0 

0 


13 


V, 


23 


0 0 

O 0 

0 0 

O 0 

0 0 


^33 

0 0 
1/4 O 
-1/4 0 

1/4 0 

-1/4 0 

1/4 0 

-1/4 0 

1/4 0 

-1/4' 0 



A 

V 


21 


/° 

0 

0 

1/J9 

-8//72 

0 

0 

0 0 

/ 

0 

1/4 

1/J9 

1/JT2 

0 

0 

0 1/4 

0 

0 

-1/4 

1/J9 

1/fn 

0 

0 

0 1/4 

0 

0 

1/4 

1/J9 

i/ffi 

0 

0 

0 - 1/4 

0 

0 

-1/4 

i/Js 

1/J72 

0 

0 

0 - 1/4 

0 

0 

1/4 

1/J9 

1/JT2 

0 

0 

0 1/4 


0 

-1/4 

1/J9 

1/Jt2 

0 

0 

0 1/4 


0 

1/4 

1/F 

I/JT 2 

0 

0 

0 

1 

1 ° 

0 

-1/4 

1//9 

1 /fn 

0 

0 

0 -1/4 



/ 0 

} 

0 

0 

0 

0 

0 

1//9 

-8//72 

0 


1 0 
i 

0 

1/4 

0 

0 

1/4 


1/,J72 

0 

1 

1 0 

0 

-1/A 

0 

0 

1/4 

1/(9 

1/(72 

0 

! 

0 

0 

-1/4 

0 

0 

-1/4 

1/J9 

1/(12 

0 

^31 = 1 

0 

0 

1/4 

0 

0 

-1/4 

Vj9 

i/fn 

0 


0 

0 

1/4 

0 

0 

1/4 

1/J9 

1/(7T 

0 


0 

0 

-1/4 

0 

0 

1/4 

1/(9 

i/.rn 

0 


1 ° 

0 

-1/4 

0 

0 

-1/4 

1/(9 

1/(72 

0 


\° 

0 

1/4 

0 

0 

-1/4 

1//9 

l/(l2 

0 

and the 

A 

submatrices V- 

L 2 ^ 

^13^ 

^^ 22 ' 

A A 

^23' ^32 

A 

and V 32 



are the same as^defined in Reference (31). 
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